function mom = compute_moments_for_calibration(calib, param, glob, options)

    %% assign parameters

    %    
    param.rho_s       = calib(1); 
    param.sigma_s     = calib(2); 
    param.sigma       = calib(3);
    param.epsilon     = calib(4)*param.sigma; 
    param.xi          = calib(5);
    param.phi_L     = calib(6);

    param.ce       = exp(-1);

    %% Set up and compute steady state
    
    [param,glob]    = setup(param,glob,options);
    
    eq_SS           = solve_eq_ss(param, glob, options);

    
    %% then compute_moments
    if eq_SS.flag == 0
      mom = compute_moments(eq_SS, param, glob, options);
    else
      mom.m1 = nan;
    end

    mom.flag = eq_SS.flag;
    
end

